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Type 1 diabetes (TID) is a multigenic disease caused by T-cell mediated destruction of the insulin producing 
pancreatic islet (3-cells. The earliest sign of islet autoimmunity in NOD mice, islet leukocytic infiltration 
or insulitis, is obvious at around 5 weeks of age. The molecular alterations that occur in T cells prior to 
insulitis and that may contribute to TID development are poorly understood. Since CD4 T-cells are essential 
to TID development, we tested the hypothesis that multiple genes/molecular pathways are altered in 
these cells prior to insulitis. We performed a genome-wide transcriptome and pathway analysis of whole, 
untreated CD4 T-cells from 2, 3, and 4 week-old NOD mice in comparison to two control strains (NOR 
and C57BL/6). We identified many differentially expressed genes in the NOD mice at each time point. 
Many of these genes (herein referred to as NOD altered genes) lie within known diabetes susceptibility 
(insulin-dependent diabetes, Idd) regions, e.g. two diabetes resistant loci, Idd27 (tripartite motif-containing 
family genes) and Iddl3 (several genes), and the CD4 T-cell diabetogenic activity locus, Idd9/ll (2 genes, 
KH domain containing, RNA binding, signal transduction associated 1 and protein tyrosine phosphatase 
4a2). The biological processes associated with these altered genes included, apoptosis/cell proliferation and 
metabolic pathways (predominant at 2 weeks); inflammation and cell signaling/activation (predominant at 

3 weeks); and innate and adaptive immune responses (predominant at 4 weeks). Pathway analysis identified 
several factors that may regulate these abnormalities; eight, common to all 3 ages (interferon regulatory 
factor 1, hepatic nuclear factor 4, alpha, transformation related protein 53, BCL2-like 1 (lies within Iddl3), 
interferon gamma, interleukin 4, interleukin 15, and prostaglandin E2); and two each, common to 2 and 

4 weeks (androgen receptor and interleukin 6); and to 3 and 4 weeks (interferon alpha and interferon 
regulatory factor 7). Others were unique to the various ages, e.g. myelocytomatosis oncogene, jun oncogene, 
and amyloid beta (A4) to 2 weeks; tumor necrosis factor, transforming growth factor, beta 1 , NFkB, ERK, and 
p38]VlAPK to 3 weeks; and interleukin 12 and signal transducer and activator of transcription 4 to 4 weeks. 
Thus, our study demonstrated that expression of many genes that lie within several Idds (e.g. Idd27, Iddl3 and 
Idd9/ ! 1 ) was altered in CD4 T-cells in the early induction phase of autoimmune diabetes and identified their 
associated molecular pathways. These data offer the opportunity to test hypotheses on the roles played by 
the altered genes/molecular pathways, to understand better the mechanisms of CD4 T-cell diabetogenesis, 
and to develop new therapeutic strategies for TID. 

® 2014 The Authors. Published by Elsevier B.V. 
This is an open access article under the CC BY-NC-ND license 
(http;//creativecommons.org/licenses/by-nc-nd/3.0/). 



1. Introduction 

Type 1 (autoimmune) diabetes (TID) is a multigenic disease that, 
in humans and NOD mice, results from T-cell mediated destruction of 
insulin producing beta cells in the pancreatic islets [1-4]. Insulitis, the 
earliest sign of autoimmune pathology in the pancreas of NOD mice, 
in our colony becomes obvious at around 5 weeks of age, at which 
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Stage accumulation of leukocytes is observed around the pancreatic 
islets. This process progressively intensifies leading to T lymphocytic 
infiltration of the pancreatic islets and eventual massive destruction 
of the insulin producing beta cells with clinical disease occurring at 
12 weeks of age or later. The molecular alterations that occur in T 
cells prior to insulitis and that may contribute to TID pathogenesis 
are poorly understood. 

It is well established that genetic predisposition is a major factor 
in the etiology of TID. The major histocompatibility (MHC) class II 
molecule H2-AS^ in NOD mice (or HLA-DQ2 and -DQ_& in humans) is 
the strongest genetic determinant for TID development [1-4]. The 
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MHC class Il-restricted CD4 T-cells are essential for the development 
of autoimmune diabetes [ 2,5 ]. To this end, ( 1 ) CD4 T-cells from spleens 
of NOD mice are reactive to pancreatic beta cell antigens [6]; (2) 
Spleen CD4 T-cells from NOD mice can transfer diabetes to young 
NOD and NOD.scid mice [7,8] ; (3) NOD mice lacking CD4 T-cells do not 
develop diabetes [9]. Reconstitution of these mice with NOD spleen 
CD4 T-cells leads to development of diabetes [5,1 0] ; and (4) transgenic 
NOD mice harboring CD4 T-cells with T cell receptor reactive to islet 
antigens develop insulitis and diabetes [8]. The chromosomal regions 
that modulate Tl D susceptibility in NOD mice are designated insulin- 
dependent diabetes regions (Idd). In conjunction with the MHC locus, 
which is required but alone insufficient for TID development, over 
20 non-MHC Idd loci also contribute to the disease process in NOD 
mice [3,4,11-14]. Although, the identities of some of the non-MHC 
Idd genes that interactively contribute to the diabetogenic process in 
CD4 T-cells of NOD mice have been revealed [5,10,15-19], most of 
these genes and/or their interactions remain unknown. 

The conventional approach in the field to understanding the 
pathogenesis of TID has been mainly via targeted analysis at the 
individual gene or loci level. Identification of all the genes that to- 
gether cause diabetes (a multigenic disease) via these approaches can 
be tortuous as each gene may only contribute weakly to the pathol- 
ogy. While these approaches have yielded useful information on how 
identified genes may interact with each other to confer disease sus- 
ceptibility and/or protection, a whole cellular and/or molecular sys- 
tems analysis (non-targeted approach) provides the opportunity to 
simultaneously interrogate the genes/pathways that are involved in 
the disease process [20]. A comprehensive understanding of these 
molecular interactions is important because it is now clear that the 
best targets for development of novel prevention and/or treatment 
interventions for complex trait diseases may not be the disease as- 
sociated genes per se but rather their interaction partners, upstream 
regulators or downstream targets, or the molecular network [21-24]. 
Thus, to gain insights into the molecular networks that might play a 
role in the diabetogenic activity of CD4 T-cells in the early induction 
phase of TID, we evaluated the transcriptomes of untreated, whole 
CD4 T-cells collected from the spleens of NOD mice in the period 
prior to overt insulitis and inferred the associated altered molecular 
networks using a suite of complementary bioinformatics tools. 

2. Materials and methods 

2. 1. Mice, sample collection and nnicroarray procedures 

Animal procedures were approved by the University of Tennessee 
(UT) Health Science Center and Veteran Affairs (VA) Medical Center 
Animal Care and Use Committee (Protocol Numbers: UT 1159/VA 
00157). Breeder mice were purchased from the Jackson Laboratory 
and housed at the VA animal facility. Spleen leukocytes were col- 
lected, as described previously 125,26] from female NOD mice at 2, 3 
and 4 weeks of age (representing the period prior to overt insulitis) 
and from two matched control strains, NOR and C57BL/6 (C57); n = 5 
for each strain and age group, except NOD 2 week, where n = 4. CD4 T- 
cells were then negatively separated by magnetic beads according to 
the manufacturer's protocol (Miltenyi Biotec). Purity was assessed by 
flow cytometry using FlTC-conjugated anti-mouse CD4 monoclonal 
antibodies (Becton Dickinson); only samples of > 90% purity were 
used in the study. Total RNA was extracted from untreated, whole 
CD4 T-cells, as described previously [26,27]. A total of 1 -1 .5 [ig of to- 
tal RNA was processed with Two-Cycle Target labeling protocol and 
hybridized on Affymetrix, Mouse430_2 expression arrays, according 
to the manufacturer's instructions. Normalization, scaling, and basic 
evaluation of the quality of the expression data from each chip were 
conducted using the GCOS software (Affymetrix), as described pre- 
viously [25,26]. The microarray data sets are available in the gene 
expression omnibus repository [GSE46600]. Microarray results were 



validated by quantitative Real-time PCR, as described previously [25] 
(Fig. SI). Expression values were normalized to glyceraldehyde-3- 
phosphate dehydrogenase (Gapdh). The target genes, primers and 
probes are listed in Table SI. 

2.2. Statistical analysis 

Statistical analysis of the microarray data was conducted as previ- 
ously described 125,26]. Filtration of the probe sets present on the chip 
array (~60,000) identified ~31,000 probe sets that had a present/ 
marginal expression flag in at least one of the samples. We then per- 
formed one-way ANOVA (at various statistical stringencies) on the 
filtered probe sets for each age separately in order to define lists of 
age-specific genes that were differentially expressed between strains. 
Lists generated ateitherp < 0.005 ("smaller" list)orp < 0.05 ("larger" 
list) adjusted with Benjamini-Hochberg multiple test correction (cor- 
responding to a false discovery rate (FDR) of 0.5% or 5%, respectively) 
were used in further analysis. Finally, we conducted hierarchical clus- 
tering on these lists as previously described 127] to identify genes that 
were uniquely differentially expressed in NOD mice relative to both 
control strains. These genes are herein referred to as NOD altered 
genes. 

2.3. Data mining analyses 

We subjected the lists of NOD altered genes to data mining using 
a suite of modern bioinformatics tools. Enriched gene ontology (GO) 
categories and KEGGs pathways were determined using WebGestalt 
Gene Set Analysis Toolkit ( http://bioinfo.vanderbilt.edu/webgestalt 
128]), as described previously 125]. To identify transcription reg- 
ulators whose binding sites were significantly over-represented 
in the promoters of the NOD altered gene sets, we used the 
PRIMA (PRomoter Integration in Microarray Analysis) program 
of the EXPANDER suite (Expression Analyzer and DisplayER; 
http : / /acgt.cs.tau.ac.il /expander/ overview.html ; Tel Aviv Univer- 
sity, Israel). PRIMA achieves this by utilizing known models for 
transcription factor binding sites. Corresponding transcription reg- 
ulators are considered to be candidate regulators of the corre- 
sponding set of genes. We used Ingenuity Pathway Analysis (IPA, 
http://www.ingenuity.com ) as described in detail previously 125] to 
identify de novo molecular networks that are associated with the NOD 
altered genes. 

3. Results 

3. 1. Genes differentially expressed in NOD CD4 T-cells at the 
preinsulitis stage of autoimmune diabetes 

We compared the whole genome mRNA expression in CD4 T-cells 
from 2-, 3-, and 4-week old NOD mice to that of two matched control 
strains, NOR and C57BL/6 (C57). NOR shares -88% of its genome 
with NOD mice, including the diabetogenic H2^^ MHC haplotype and 
several important non-MHC TID susceptibility loci 11-4,29]. C57, on 
the other hand, is a more genetically distantly related strain to NOD. 
Yet, like C57, NOR is both insulitis- and diabetes-free. 

We identified 362, 982, and 581 probe sets (genes) with highly 
significant expression differences between strains (p < 0.005, 
Benjamini-Hochberg; FDR of 0.5%) at 2, 3, and 4 weeks, respectively. 
As expected, the majority of these genes had a similar pattern of ex- 
pression in NOD and NOR compared to C57 (Fig. 1 ). The focus of our 
study was to identify genes in NOR mice whose expression was similar 
to that in C57 but different from that in NOD, i.e. genes differentially 
expressed in NOD relative to both NOR and C57, and herein referred to 
as NOD altered genes. These constitute prospective candidate genes 
for protection of NOR mice against diabetes. Thus, we identified a total 
of 58, 115, and 65 probe sets whose expression was altered in NOD 
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at 2, 3 and 4 weeks, respectively, compared to both controls (clus- 
ters of lower or higher expression in NOD are indicated by arrows in 
Fig. 1 ). These represented 56, 1 07, and 60 different genes, respectively, 
shown in Tables 1 -4. The great majority of these genes were of lower 
expression in NOD mice compared to controls: ~70% at 2 weeks; ~72% 
at 3 weeks; and ~87% at 4 weeks. Twenty-five genes (72% of which 
was of lower expression in NOD mice) were common to all 3 ages. 
Results for qRT-PCR validation of eight genes are shown in Fig. SI (3 
of higher expression in NOD mice: protein tyrosine phosphatase 4a2 
(Ptp4a2), biogenesis of organelles complex-1 , subunit 6 (Blocl s6, pal- 
lidin) and transmembrane protein 87A (Tmem87a); and 5 of lower 
expression in NOD mice: tripartite motif-containing 5 (TrimS), tri- 
partite motif-containing 1 2A (Triml 2a), Gatm (glycine amidinotrans- 
ferase (L-arginine:glycine amidinotransferase)), lymphocyte antigen 
6 complex, locus CI (Ly6cl ), and receptor transporter protein 4 (Rtp4) 
compared to control mice). The major selected functional categories 
(as determined by the authors based on information obtained from 
literature searches) included immune response, apoptosis/cell pro- 
liferation, transcription, zinc-ion binding, protein/nucleic acid bind- 
ing, "Other enzymes" for enzymes not already assigned to a specific 
category, and "Other" for genes that did not fall in one major cate- 
gory. Probe sets of unknown identification/function are listed under 
"Unknown". Based on these categories, a higher percentage of NOD 
altered genes (~22%) were involved in apoptosis/cellular prolifera- 
tion at each of 2 or 3 weeks than at 4 weeks (12%). Conversely, a 
higher percentage of genes (22%) were involved in immune response 
at 4 weeks than at 2 or 3 weeks (4.1% and 12%, respectively). Overall, 
a high percentage of genes at all 3 ages were involved in enzymatic 
activity, including several kinases/phosphatases at 3 weeks. More- 
over, a large number of the CD4 T-cell NOD altered genes lie within 
known Idds, 52%, 47% and 38% at 2, 3, and 4 weeks, respectively 
(Tables 1-4). A total of 14 (25%), 20 (-19%) and 10 (-17%) genes 
at 2, 3, and 4 weeks, respectively, lie within lddl3, a region located 
on chromosome (Chr) 2 and previously identified to confer resis- 
tance to diabetes in NOR mice [4,15-17,30]. These genes included 
6 that were common to all 3 ages: Blocl s6, Trp53bpl (transforma- 
tion related protein 53 binding protein 1), Tmem87a, Ctdspl2 (CTD 
(carboxy-terminal domain, RNA polymerase II, polypeptide A) small 
phosphatase like 2), Gatm and Raly (hnRNP-associated with lethal 
yellow). Two genes (Khdrbsl - KH domain containing, RNA binding, 
signal transduction associated 1 and Ptp4a2, both altered at 2 and 

3 weeks, but only Khdrbsl at 4 weeks) lie within ldd9/l I, a region on 
Chr4 that also confers resistance to diabetes in NOR. This region has 
also been demonstrated previously by several groups to regulate the 
diabetogenic activity of CD4 T-cells [5,1 0,1 8,1 9]. 

3.2. Gene ontology categories and KEGG pathways of the CD4 T-cell 
NOD altered genes 

To gain further insights into the biological processes associated 
with the NOD altered genes, we performed Gene ontology (GO) and 
KEGG pathway analyses. Because GO analysis is better suited for larger 
gene lists, we also analyzed gene lists generated at a slightly lower 
statistical stringency. Analysis atp < 0.05, Benjamini-Hochberg (FDR 
of 5%) identified 134, 252, and 185 NOD altered genes at 2, 3, and 

4 weeks, respectively (Tables S2, S3 and S4, respectively). The topmost 
ranked major biological process of these gene lists was metabolism 
(Table 5). Although common to all 3 ages, this functional category 
was most significantly enriched at 2 weeks. Several biological pro- 
cesses unique to 3 weeks were also identified, including biopolymer 
modification, localization and transport, and T cell activation. Un- 
der "molecular function", hydrolase activity and binding were the 
common top ranked categories. "Intracellular" was the most signif- 
icantly enriched cellular component for all the 3 lists. Nucleus and 
cytoplasm/endoplasmic reticulum were uniquely enriched at 2 and 
3 weeks, respectively. The immune system process ranked at the 



top of the biological process at 4 weeks in the analysis of the gene 
list generated at the higher stringency level (p < 0.005, Benjamini- 
Hochberg; FDR of 0.5). Findings of the GO analyses were consis- 
tent with those from our "manual" literature searches. We then per- 
formed KEGG pathway analysis to identify well characterized molec- 
ular pathways that were significantly over-represented in the gene 
lists of NOD altered genes. Results of the smaller and larger gene lists 
of NOD altered genes were similar; only those of the smaller lists 
are presented (Table 6). Consistent with the GO analysis, the pre- 
dominantly enriched category was metabolic pathways, which still 
was most highly significantly enriched at 2 weeks as compared to 
the other two ages. Five NOD altered genes common to all 3 ages 
were identified in these metabolic pathways: Enpp3 (ectonucleotide 
pyrophosphatase/phosphodiesterase 3), Ndufs5 (NADH dehydroge- 
nase (ubiquinone) Fe-S protein 5), GalntlO (UDP-N-acetyl-alpha-D- 
galactosamine:polypeptide N-acetylgalactosaminyltransferase 1 0), 
Prim2 (DNA primase, p58 subunit) and Gatm. These findings suggest 
that these genes may contribute to the metabolic abnormalities that 
affect the immune system and predispose NOD mice to autoimmune 
diabetes. Overall, these data suggest that CD4 T-cells from NOD mice 
already have a defect in metabolism (most prominent at 2 weeks), 
cellular activation and endoplasmic reticulum function (both evident 
at 3 weeks) and T-cell/immune response (evident at 3-4 weeks) at 
the preinsulitis stage. 

3.3. Transcriptional regulatory pathways of the CD4 T-cell NOD altered 

genes 

The 4 topmost significantly enriched transcription factor bindings 
sites (Table 7) were for androgen receptor (Ar), significantly enriched 
at 2 and 4 weeks; Interferon regulatory factor 1 (Irfl ), significantly en- 
riched at all 3 ages; and Interferon regulatory factor 7 (Irf7) and Inter- 
feron sensitive response elements (ISRE) both significantly enriched 
at 3 and 4 weeks. ISRE are present in the promoters of interferon 
stimulated genes (ISGs), also known as antiviral or innate immune 
response genes. All 4 binding sites were most highly significantly en- 
riched at 4 weeks. Interestingly, each putative transcription regulator 
correlated with a different set of NOD altered genes at the various 
stages it was significantly enriched, including both age-common and 
age-specific NOD altered genes. This suggests a dynamic (possibly co- 
ordinated) regulation of gene expression. The promoter analysis data 
suggest that expression of the NOD altered genes was significantly 
more likely to be regulated by Ar, Irfl, lrf7, and type I interferon than 
by other transcription regulators. 

3.4. Ingenuity Pathway Networks of the CD4 T-cell NOD altered genes 

In addition to the transcriptional regulatory pathway analysis, we 
performed Ingenuity pathway analysis (IPA; 17,18) to gain further 
insights into the genes/molecules that may play a role in regulating 
the expression of the NOD CD4 T-cell altered genes and/or molec- 
ular pathways. IPA analyses of the 3 "smaller" lists of NOD altered 
genes (Tables 1-4) generated several networks each. The major bio- 
logical functions significantly represented by the networks included 
cell cycle, cellular growth/proliferation, and cell death at 2 weeks; 
molecular transport, cell-to-cell signaling and interaction, and im- 
mune system development and function at 3 weeks; and antiviral 
(innate immune) function and immune response at 4 weeks. The 
merged networks for each dataset clustered around several central 
genes (Figs. 2-4). To facilitate an objective comparison between net- 
works and to gain further insights, we ranked the central genes ac- 
cording to the total number of connections linked to each one of 
them - the more the number of connections, the higher the rank. 
The top 15 central genes and their directly linked focus genes are 
shown in Table 8. This analysis revealed that seven central genes 
were common to all 3 ages, including transcription factors, hepatic 
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Table 1 

Genes differentially expressed in CD4 T-cells from 2 week-old NOD mice. 
Probe set ID Gene symbol Fold change Adjusted p-value Chromosome (Chr) TID susceptibility region 



Immune response 



1418642 at (L) 


Lcp2 


-9.65 


0.00048 


Chrll 


- 


1457088_at (H) 


Pldn 


5.48 


0.00128 


Chr2 


/ddl3 


Apoptosis/cell proliferation 












1440493 at (L) 


GalntlO 


-37.09 


0.000217 


Chrll 


- 


1457812_at (L) 


Trp53bpl 


-30.33 


0.000291 


Chr2 


/ddl3 


1438462_x_at (L) 


Khdrbsl 


-7.42 


0.000662 


Chr4 


Idd9/ll 


141 7714 j(_at(H) 


Hba-al 


5.36 


0.0031 


Chrll 


- 


1449716_s_at (H) 


Nrdl 


2.77 


0.00321 


Chr4 


- 


1439650_at(L) 


Rtn4 


-5.74 


0.00296 


Chrll 


- 


1440358_at(H) 


ArhgeflS 


4.82 


0.00184 


Chrll 


- 


mziSjo 1 _x_dU {Lj 


(Jill t1 


yl 1 1 


U.UU lU'i 


Llll 1 1 




1444890_at(L) 


Mprip 


-4.08 


0.00368 


Chrll 




1449052_a_at(L) 


Dnmt3b 


-3.24 


0.00048 


Chr2 


Iddl3 


1433681 _x_at(L) 


Capn3 


-2.30 


0.00214 


Chr2 


Iddl3 


Transcription 












1458094 at (L) 


Zfp407 


-6.66 


0.00224 


Chrl8 




1458274_at (H) 


Zfp69 


3.06 


0.00288 


Chr4 


- 


1446147 _at(H) 


Rbm39 


4.15 


0.00129 


Chr2 


Md!3 


Transport 












1455735_at(H) 


Apls3 




U.UUUj iz 




IQClD.HU/ D.H 


1417963_at(H) 


PItp 


3.00 


0.00419 


Chr2 




1420897_at(L) 


Snap23 


-1.81 


0.000668 


Chr2 


lddJ3 


Zinc-ion binding 












1437432_a_at (L) 


Triml2a 


-789.41 


2.45E-08 


Chr7 


Idd27 


1443858_at (L) 


Trim/12c (Trim5) 


-238.91 


3.05E-07 


Chr7 


Idd27 


Protein/nucleic acid binding 










1424454 at (H) 


Tmem87a 


13.40 


0.00114 


Chr2 


IddU 


1455863 _at (L) 


Spata511 


-5.40 


0.00217 


Chr2 


- 


1419276_at(H) 


Enppl 


2.07 


0.00105 


ChrlO 


- 


1453065_at(H) 


Aldh5al 


1.78 


0.00337 


Chrl3 


- 


Otiier enzymes 












1442424_at (L) 


Ctdspl2 


-46.26 


0.000831 


Chr2 


IddU 


1418035 a at (L) 


Prim2 


-11.25 


5.55E-05 


Chrl 


Idd26 


1423569_at (L) 


Gatm 


-7.41 


n nnni 1 1 
U.UUU IIZ 


Chr2 


Idd13 


1453009 _at (L) 


Cpm 


-6.18 


0.00158 


ChrlO 




1416494 at (L) 


Ndufs5 


-5.21 


0.000538 


Chr4 


- 


1427943_at (H) 


Acyp2 


2.75 


4.07E-05 


Chrll 


- 


1427302_at (H) 


Enpp3 


2.66 


0.0027 


ChrlO 


- 


1417826_at(L) 


Akrlel 


-6.52 


0.000273 


Chrl 3 


- 


1455219_at(L) 


Slxlb 


-5.00 


0.00184 


Chr7 


Not Assignei 


1442466_a_at(L) 


Ppip5kl 


-414 


0.00211 


Chr2 


Mdl3 


1444377_at(L) 


Psmb2 


-3.41 


0.0028 


Chr4 


Mdll 


1435129_at(H) 


Ptp4a2 


2.37 


0.00487 


Chr4 


;dd9/!I 


1454772_at(L) 


Snrnp200 


-2.32 


0.00189 


Chr2 


MdI3 


1415878_at(L) 


Rrml 


-1.90 


0.00407 


Chr7 


- 


1451998_at(L) 


Taspl 


-1.78 


0.000499 


Chr2 


idd 13 


Otiier (unknown transcripts are indicated below in the legend) 








1444741_at(L) 


Dock2 


-144.49 


3.18E-06 


Chrll 


- 


1436061 at (L) 


Chafla 


-21.85 


0.000538 


Chrl 7 




1442824_at (L) 


Raly 


-13.69 


0.000313 


Chr2 


Idd13 


1452426_x_at (L) 


LOC433762 


-12.00 


0.000549 


Chr4 




1452359 at (L) 


Rein 


-2.53 


0.000662 


Chr5 




1457822_at(L) 


Tmeml31 


-13.55 


0.00438 


Chrl 


Idd26 


1445214_at(L) 


Rex2 


-9.26 


0.000783 


Chr4 


ldd92 


1456635_at(L) 


SpllO 


-8.01 


0.00156 


Chrl 


Idd5.4a/5.4 


1428587_at(L) 


Tmem41 b 


-436 


0.000289 


Chr7 


ldd27 


1424509 _at(H) 


Cdl77 


4.03 


0.00284 


Chr7 




1435792_at(H) 


Csprs 


3.86 


0.00199 


Chrl 


ldd5.4a/5.4 


1458684_at(H) 


Ssl8 


2.77 


0.00407 


Chris 


Idd21.3 


1424721 _at(L) 


Mfap3 


-2.54 


0.000895 


Chrll 




1425331 _at(L) 


Zfpl06 


-2.10 


0.00183 


Chr2 


Idd13 



Statistical analysis (1-way-ANOVA; p < 0.005, Benjamini-Hochberg) followed by hierarchical clustering identified 58 differentially expressed probe sets (representing 56 different 
genes) at 2 weeks; L and H, respectively, indicate lower or higher expression in NOD mice compared to control mice, NOR and C57BL/6. Two unknown transcripts (both on 
Chr4) included: 9930104L06Rik (L) and 2510039O18Rik (L; Idd92). Fold change (FC) was calculated by ratio of means of expression in NOD mice versus controls. Dashes indicate 
those genes are not located within a known TID susceptibility region (/dd); all indicated Idds (except one, Idd5.4a/5.4) were identified as conferring resistance to diabetes 
( http://www.tl dbase.org ). Genes highlighted in bold font were differentially expressed at all 3 ages, 2, 3 and 4 weeks. 



34 



D.N. Kakoola et al. / Results in Immunology 4 (2014) 30-45 



Table 2 

Genes differentially expressed in CD4 T-cells from 3 week-old NOD mice (continued in Table 3). 





Gghg symbol 


Fold change 


Adjusted p-value 


Chromosome (Chr) 


TID susceptibiHty region 


Immune response 












141sb4Z_at (L) 


Lcp2 


-12.20 


U.UUJJ9 


Chrll 




1457Uoo_at (H ) 


Pldn 


5.67 


n i\ni\Aot 
U.UUU4ol 


Chr2 


lddl3 


1 A AO !zrzr\ 

144oDDU_at [L) 


Lbp 


-3,98 


U.UUUl /b 


Chr2 


lddl3 


14JbDzy_at [L) 


Ifitllb 


-2.91 


A A AT A T 

U.UUz4z 


Chrl9 




14DUJy4_a_dC l^n J 


Inppll 


2.76 


A AAQ 1 Q 
U.UU J 1 O 


Chr7 




14DU/oJ_at (LJ 


Ifitl 


-2.73 


A AAOAC 

u.uujyb 


Chrl9 




14z4// b_at (LJ 


Oasla 


—2.20 


0.001 15 


Chr5 




1 /1 1 O C on -.^ / T \ 

141oDoU_at (LJ 


Rtp4 


—2.10 


A AAO 1 "7 
U.UU31 / 


Chrl6 




144oy4U_at [L) 


Trim21 


-2.03 


A AAOAA 

u.uujyy 


Chr7 


!dd27 


14J jjfciu_ar [L) 


Itgal 


—1.84 


U.UUz4z 


Chr7 


Not Assigned 


i4DzyDb_a_at [Lj 


Ifi27 


-1.74 


A AA/I O A 


Chr6 




• 

Apoptosis/celi proliferation 












145781Z_at (Lj 


Trp53bpl 


-16.30 


i\ i\t\i\r\tiA 
0.000^64 


Chr2 


Iddl3 


■tAAnAtfy -t*- fi \ 

I4404yj_at (Lj 


GalntIO 


-11.41 


U.UU074o 


Chrll 




1417714_X_at (rlj 


Hba-al 


10.69 


i\ r\c\r\-t cc 
U.UUUl 5b 


Chrll 




14Jo4fa2_X_at (LJ 


KtiarDsl 


-5.59 


\i.\i\3\iZ\ 1 


Chr4 


ma^f 1 1 


l44y/lD_s_at [H) 


Nrdl 


2.36 


0.00133 


Chr4 




14Dzb/ /_at (n J 


Pnptl 


13.80 


0.0012 


Chrll 




144UJ 1 y_ar [Lj 


Mef2a 


-13.53 


U.UUzz / 


Chr7 




1444oyU_at (LJ 


Mprip 


-9.37 


A r\c\A on 
U.UU43y 


Chrll 




1439650_at(L) 


Rtn4 


-5.24 


0.00318 


Chrll 


I 


1433681 _x_at(L) 


Capn3 


-2.89 


0.000442 


Chr2 


lddl3 


1 ylOQQI Q it- (\ \ 

mztso 1 y_ar [Lj 


Maprcl 


*> 70 

—£.,IO 


U.UUZ4Z 


Lnrz 


IQQU 


i44yujz_a_ac i^lj 


unrntJD 


—2.77 


1.06E-05 


Chr2 


Iddl3 


143bb4/_at (hi J 


Ttbl<2 


2.69 


A AA1C A 

U.UUzb4 


Chr2 


!ddl3 


i4iyzby_at (h j 


Dut 


2.54 


0.0016 


Chr2 


!ddl3 


1444DUU_at (LJ 


Ahsal 


—2.48 


0.0042 


Chrl2 




14JbUzb_at (n J 


CcdcSSa 


2.40 


A AAACA"? 

u.uuuby / 


Chrll 




I44iyj/_s_at [Lj 


Pinkl 


-2.11 


A AA1 CA 

u.uuiby 


Chr4 




144olo4_at (LJ 


rkDpla 


—1.94 


A AA/I A A 

U.UlJ44y 


Chr2 


Idd'!3 


14zzoUo_S_at (LJ 


Dock2 


—1.84 


A AAAAAA 

u.uuuyyy 


Chrll 




14d /ol3_at (n J 


IrpbJDpl 


1.62 


A A A /I A 

U.UU44 


Chr2 


!ddl3 


144oUz/_at (LJ 


Ncoa3 


—1.47 


A AAT7"7 

U.UUz / / 


Chr2 




1419562_at(L) 


Birc6 


-131 


0.00393 


Chrl7 


_ 


Tra nscription 












1458l)94_at (L) 


Zfp407 


-8.70 


0.0uUs7S 


Chris 


Idd21.1 


1458Z74_at (HJ 


Zfp69 


3.25 


U.UUzb4 


Chr4 




i44yjyz_ar (n j 


Tcfl 5 


10.64 


n nnnoQQ 
U.UUUyiSo 


Chr2 


!ddl3 


i4Dyuzb_at (LJ 


Snwl 


-10.00 


A AA-| CI 

U.UUlbz 


Chrl2 




141 7961 _a_at(L) 


Trim30 


-7.39 


2.84E-05 


Chr7 


ldd27 


1445214_at(L) 


Zfp715 


-4.30 


0.000612 


Chr7 


- 


1^4/ / uj_x_at ^Lj 






J.4/t-Uj 


Lnr^ 




144z3Db_at (LJ 


iVIax 


—2.24 


U.UU4yi 


Chrl2 




14zb/bb_at (LJ 


Commd? 


-2.09 


r> Ann /I 
U.UUU4bb 


Chr2 


lddl3 


143641 6_x_at(L) 


Fxcl 


-1.81 


6.23E-05 


Chr7 


ldd27 


1450350_a_at(L) 


Jdp2 


-1.39 


0.00475 


Chrl2 




Transporters 












146061 7_s_at(L) 


Rab6b 


-15.69 


0.000612 


Chr9 




1458426_at(LJ 


Kiflb 


-3.72 


0.000704 


Chr4 


ldd9.2 


1455735_at(H) 


Apls3 


3.09 


0.00386 


Chrl 


ldd5.4a/5.4 


142461 5_at(LJ 


Pgap2 


-2.18 


0.0019 


Chr7 


ldd27 


1424211_at(H) 


Slc25a33 


1.83 


0.0012 


Chr4 




Zinc-ion binding 












1437432_a_at (L) 


Triml2a 


-252.63 


1.20E-08 


Chr7 


ldd27 


1443858_at (L) 


Triml2c (TrimS) 


-140.08 


5.39E-07 


Chr7 


Idd27 


1435665_at(LJ 


Trim30d 


-3.16 


0.00227 


Chr7 


ldd27 



Statistical analysis (1-way-ANOVA; p < 0.005, Benjamini-HochbergJ followed by hierarchical clustering identified 115 differentially expressed probe sets (representing 107 
different genes) at 3 weeks; L and H, respectively, indicate lower or higher expression in NOD mice compared to control mice, NOR and C57BL/6. This table shows a partial list of the 
genes; the remaining genes are shown in Table 3. Fold change (PC) was calculated by ratio of means of expression in NOD mice versus controls. Dashes indicate those genes are not 
located within a known TID susceptibility region (Idd); all indicated /dds (except one, ldd5.4a/5.4) were identified as conferring resistance to diabetes ( http;//www.tldbase.org ). 
Genes highlighted in bold font were differentially expressed at all 3 ages, 2, 3 and 4 weeks. 
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I NOD I NOR IC57BL/6 B I NOD I NOR I C57BL/6 0 | NOD | NOR | C57BL/6 




Fig. 1. Hierarchical clusters of genes wliose expression was altered in CD4 T-cells. 362 (A). 982 (B) and 581 (C) genes were differentially expressed between strains at 2, 3, and 
4 weeks of age, respectively. The lists were identified by a one-way ANOVA of ~31,000 filtered probe sets at p < 0.005, with Benjamini-Hochberg multiple test correction. A total 
of 58, 115, and 65 probe sets were differentially expressed in NOD relative to both controls (NOR and C57BL/6) at 2-, 3- and 4-weeks, respectively; clusters of these NOD altered 
genes are indicated by arrows. The color intensity of the rectangles representing each gene for each sample (n = 5 for each strain/age, except NOD 2 week, where n = 4) indicates 
the degree of increase (red) or decrease (blue) of the gene expression signal relative to the mean signal intensity (yellow). 



nuclear factor 4, alpha (Hnf4a) and transformation related protein 53 
(Trp53); anti-apoptotic factor, BCL2-lii<e 1 (Bcl2n); cytokines, inter- 
feron gamma (Ifng), interleukin 4 (114) and interleukin 15 (1L15); and 
the pro-inflammatory factor, prostaglandin E2. The pro-inflammatory 
cytokine, interleukin 6 (116) was common to ages 2 and 4 weeks. Other 
genes were unique to the respective age groups, e.g. to 2 weeks: tran- 
scription factors, myelocytomatosis oncogene (Myc), Jun oncogene 
(Jun), and amyloid beta (A4) precursor protein (App); to 3 weeks; 
the pro-inflammatory cytokine, tumor necrosis factor (Tnf), the anti- 
inflammatory cytokine, Tgfbl (transforming growth factor, beta 1 ), 
and cell signaling molecules NFkB, ERK, p38MAPl<, and insulin-like 
growth factor 1 (Igfl ); and to 4 weeks: cytokines, interferon alpha and 
1L12 (interleukin 12), the transcription factor signal transducer and 
activator of transcription 4 (Stat4), and Rous sarcoma oncogene (Src). 
Interestingly, the central genes Bcl211 and Src are also located within 
Iddl3. The biological processes associated with the central genes that 
were common to all 3 ages included apoptosis/cellular proliferation, 
Thl-Th2 balance, cytokine signaling, and inflammation. Those as- 
sociated with the age-specific central genes included apoptosis/cell 
proliferation (at 2 weeks), cell signaling/cellular activation and in- 
flammation (at 3 weeks), and innate immunity/link to adaptive im- 
mune response (at 4 weeks). Similar to the promoter analysis, each of 
the common central genes was connected to a different set of focus 
genes at the various ages, consisting of both age-common and age- 
specific genes, again suggesting a dynamic coordinated regulation of 
the molecular processes over time. These data suggest that abnor- 
malities in apoptosis/cellular proliferation predominate at 2 weeks 
of age, while those in cell signaling/cellular activation and inflam- 
mation predominate at 3 weeks of age. Abnormalities in the innate 
immune response and its link to adaptive immunity predominate at 
4 weeks of age. 



4. Discussion 

Type 1 diabetes is a multigenic disease whose study would benefit 
from a non-targeted approach to simultaneously identify the genes/ 
pathways that are involved in the disease process [20]. Here, we em- 
ployed a novel 3-way analysis where NOD was compared to two 
diabetes-resistant strains, C57BL/6 and NOR. Comparison with NOR 
mice (a strain that is ~88% identical to NOD mice, including shar- 
ing the diabetogenic MHC haplotype) afforded us the ability 
to identify gene expression changes (and their associated molecular 
networks) that might contribute to TID susceptibility in NOD or, con- 
versely, resistance in NOR mice in the context of a permissive MHC 
haplotype. This study confirmed as well as advanced our recently 
published studies on unfractionated spleen leukocytes [25] in that, 
in addition to using the NOR control strain, we focused on a specific 
leukocyte subset (CD4 T-cells) and included an additional time point 
(3 weeks of age), all of which allowed us to identify novel altered 
pathways. 

NOR mice remain diabetes-free due to resistance alleles within the 
~12% portion of their genome derived from BKs [4,17,29,30]. These 
BKs-derived genomic regions are present on chromosomes 1 , 2, 4, 5, 7, 
10, 11, 12, 14 and 18 [17,29,30]. Unsurprisingly, virtually all the NOD 
altered genes identified in the current study are located on these chro- 
mosomes (Tables 1-4). Strikingly, a large number of these genes lie 
within known diabetes susceptibility regions (Tables 1-4). Linkage/ 
congenic studies have provided evidence for the presence of NOR re- 
sistance genes on Chrl {Idd5.2), Chr2 {Iddl3) and Chr4 {IddQ/n). and 
potentially on Chrll [15-17,30]. However, several studies also sug- 
gest that genes in other genetic regions distinguishing NOR from NOD 
may also act interactively with genes in these loci to protect NOR mice 
from autoimmune diabetes [5,10,30], thus providing further support 
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Tables 

Genes differentially expressed in CD4 T-cells from 3 weel<-old NOD mice (continued from Table 2). 



Probe set ID 



Gene symbol 



Fold change 



Adjusted p-value 



Chromosome (Chr) 



TID susceptibility region 



Protein/nucleic acid binding 



1424454 
1455863. 

1459957. 
1429247. 
1458684. 
1419276. 
1416559. 
1440416. 
1429337. 
1428308. 



at(H) 
.at(L) 

.at(H) 
.at(L) 
at(H) 
.at(H) 
.at(L) 
at(H) 
.at(L) 
.at(L) 



Tmem87a 
SpataSll 

Tnrc6a 

Anxa6 

Ssl8 

Enppl 

Rrp8 

Usp46 

Tmem87b 

Pdrgl 



13^1 
-10.62 

5.53 

-4.00 

2.56 

2.48 

-232 

2.07 

-1.52 

-1.36 



0.00264 
7.71E-06 

0.000595 

0.000999 

0.00101 

0.000575 

0.00226 

0.00153 

0.00159 

0.00361 



Chr2 
Chr2 

Chr7 

Chrll 

Chris 

ChrlO 

Chr7 

Chr5 

Chr2 

Chr2 



/ddl3 

Not Assigned 

Md4.3 

ldd21.3 

ldd27 

Iddl3 
lddl3 



Otiier enzymes 
1442424_at (L) 
1453009_at (L) 
1427302 at (H) 
1418035_a_at (L) 
1423569_at (L) 
1416494 at (L) 
1427943_at (H) 
1444377_at(L) 
1449862_a_at (H) 
1451277_at(L) 
1435129_at(H) 
142061 3_at(L) 
1455075_at(L) 
1419125_at(L) 
1453343 _s_at (H) 



Ctdspl2 

Cpm 

Enpp3 

Prim2 

Catm 

Ndufs5 

Acyp2 

Psmb2 

Pi4k2b 

Zadh2 

Ptp4a2 

Ptp4a2 

Pigv 

Ptpnl8 

Vrk2 



-53.79 

-16.36 

3.45 

-10.26 

-8.19 

-6.65 

3.91 

-2.86 

2.85 

-2.50 

2.47 

-2.22 

-2.15 

-1.99 

1.97 



0.000129 

0.000671 

0.0016 

l.llE-05 

7.47E-05 

1.44E-05 

0.000117 

0.00325 

0.00324 

0.000108 

0.0033 

0.000306 

0.00133 

0.00395 

0.00454 



Chr2 

ChrlO 

ChrlO 

Chrl 

Chr2 

Chr4 

Chrll 

Chr4 

Chr5 

Chris 

Chr4 

Chr4 

Chr4 

Chrl 

Chrll 



Iddl3 



Idd26 
Iddl3 



Iddll 

ldd21.1 
ldd9/n 
Md9/n 

ldd26 



Other 

1444741 _at(L) 

1436061 at (L) 

1442824_at (L) 

1452426_x_at (L) 

1452359 at (L) 

1444878_at(L) 

1428318_at(H) 

1424133_at(H) 

1432391 _at(L) 

14341 12_at(L) 

1428587_at(L) 

1456685_at(L) 

1451477_at(H) 

1444714_at(L) 

1453465_x_at (L) 

1431259_at(H) 

1424721_at(L) 

1445186_at(H) 

1416935_at(L) 

Unloiown 

1428604_at(H) 
1434550_at(L) 
1436388_a_at (L) 
1440214_at(L) 
1429203_at(H) 
1438646_x_at (L) 
145841 9_at(L) 



Dock2 
Chafla 
Raly 

LOC433762 
Rein 

Nfl 

Smgc 

Tmem9S 

Ccdc21 

Lphn2 

Tmem41b 

HMP19/Nsg2 

Znf41-ps 

Dcdc2b 

Ppplrl4c 

Adal 

Mfap3 

Stc2 

Trpv2 



2610305D13Ril< 

3830406C13Rik 

3830406C13Rik 

A630001G21Rik 

241 007612 IRik 

2510039O18Rik 

E130215H24Rik 



-46.41 

-24.94 

-23.88 

-22.39 

-2.84 

-10.91 

7.02 

7.32 

-5.10 

-3.65 

-3.57 

-3.58 

2.90 

-2.77 

-2.30 

2.28 

-1.93 

1.78 

-1.61 



28.19 

-16.63 

-10.85 

-8.76 

3.10 

-2.69 

-2.11 



2.48E-05 

0.000216 

0.000306 

0.00114 

0.00152 

0.000999 

0.0019 

4.09E-05 

0.000648 

0.00201 

1.06E-05 

0.000655 

2.66E-05 

0.00166 

0.00159 

0.0044 

0.00143 

0.00367 

0.00289 



0.00125 

0.00157 

0.0025 

0.0044 

5.59E-05 

0.00214 

0.00254 



Chrll 

Chrl 7 

Chr2 

Chr4 

Chr5 

Chrll 

Chrl 5 

Chrll 

Chr4 

Chr3 

Chr7 

Chrll 

Chr4 

Chr4 

ChrlO 

Chr2 

Chrll 

Chrll 

Chrll 



Chr4 

Chrl4 

Chrl4 

Chrl 

Chr9 

Chr4 

Chr2 



IddU 



ldd27 
ldd9.2 

lddl3 



ldd9.2 



ldd5.4a/5.4 
ldd2 
ldd9.2 
Md!3 



Statistical analysis (1-way-ANOVA; p < 0.005, Benjamini-Hochberg) followed by hierarchical clustering identified 115 differentially expressed probe sets (representing 107 
different genes) at 3 weeks; L and H, respectively, indicate lower or higher expression in NOD mice compared to control mice, NOR and C57BL/6. This table shows a continuation of 
the gene list from Table 2. Fold change (FC) was calculated by ratio of means of expression in NOD mice versus controls. Dashes indicate those genes are not located within a known 
TID susceptibility region (/dd); all indicated /dds (except one, ldd5.4a/5.4) were identified as conferring resistance to diabetes ( http://www.tldbase.org ). Genes highlighted in 
bold font were differentially expressed at all 3 ages, 2, 3 and 4 weeks. 



for the significance of investigations at the whole cellular or molec- 
ular systems level. In further support for whole molecular systems 
studies, F2 segregation studies and subsequent subcongenic analy- 
ses have found that lddl3 (Chr2) and Idd9/ll (Chr4) each consist of 
multiple genes that contribute to NOR resistance [10,16,17,30]. Our 
study revealed that a total of 26 different NOD CD4 T-cell altered 
genes lie within lddl3, including 6 (Blocls6, Trp53bpl, Tmem87a, 
Ctdspl2, Gatm and Raly) that were common to the 3 ages studied 



(Tables 1-4). Interestingly, two central genes in the IPA networks, 
Bcl211 and Src (Table 8), also lie within this locus. Bcl211 is linked to 3 
genes, GalntlO, Rtn4 (reticulon 4) and Pnptl (polyribonucleotide nu- 
cleotidyltransferase 1) that are located on Chrll, while Src is linked 
to 2 genes, Paqr7 (progestin and adipoQ receptor family member VII) 
and Khdrbsl that are located on Chr4. Of note, Khdrbsl is located 
within a known Idd9/ll (discussed below) while the 3 Bcl211 linked 
genes and Paqr7 are not located within a known Idd, and thus may 
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Table 4 

Genes differentially expressed in CD4 T-cells from 4 week-old NOD mice. 
Probe set ID Gene symbol Fold change Adjusted p-value Chromosome (Chr) TID susceptibility region 



Iiiiiikune response 












1418642 at (L) 


Lcp2 


-20.53 


7.58E-06 


Chrll 


_ 


1457088_at (H) 


Pldn 


4.60 


0.0027 


Chr2 


/ddl3 


1421571_a_at(L) 


Ly6cl 


-9.30 


8.63E-07 


Chrl5 


- 


1435529_at(L) 


Ifitl 


-3.69 


9.11E-06 


Chrl9 


- 


1419599 _s_at(L) 


Ms4a6d 


-3.46 


1.85E-05 


Chrl9 


- 


1436058_at(L) 


Rsad2 


-2.913 


0.00244 


Chrl2 


- 


1450783 _at(L) 


ffltU 


-2.58 


0.00101 


Chrl9 


- 


1449025_at(L) 


Ifit3 


-2.39 


0.00306 


Chrl9 


- 


1424775_at(L) 


Oasla 


-2.33 


0.0015 


Chr5 


- 


1422903_at(L) 


Ly86 


-2.26 


0.00158 


Chrl3 


Iddl4 


1418580_at(L) 


Rtp4 


-2.23 


0.000429 


Chrl6 


- 


Apoptosis/cell proliferation 












1457812_at (L) 


Trp53bpl 


-20.92 


0.00426 


Chr2 


Idd13 




Udint 1 u 






unri 1 




1438462j(_at (L) 


Khdrbsl 


-9.83 


2.51E-06 


Chr4 




1417714_x_at(H) 


Hba-al 


5.67 


0.00182 


Chrll 


- 


1449716ji at (H) 


Nrdl 


2.62 


1.50E-05 


Chr4 


- 


1428371 _at(H) 


Ttbk2 


1.99 


0.0025 


Chr2 


Iddl3 


Transcription 












1458094 at (L) 


Zfp407 


-5.59 


0.0015 


Chrl8 


Idd21.1 


1458274_at (H) 


Zfp69 


2.97 


0.00361 


Chr4 


- 


145902S_at(L) 


Snwl 


-5.40 


6.75E-05 


Chrl2 


- 


1447703_x_at(L) 


Zfp593 


-4.20 


0.00209 


Chr4 


- 


1417961_a_at(L) 


Trim30 


-3.58 


0.000703 


Chr7 


Idd27 


1421519_a_at(L) 


Zfpl20 


-2.30 


0.000974 


Chr2 


lddl3 


Transporters 












1434914_at(L) 


Rab6b 


-36.11 


0.000308 


Chr9 


_ 


1458426_at(L) 


Kiflb 


-3.49 


0.00345 


Chr4 


Idd9.2 


1421098_at(L) 


Stapl 


-1.39 


0.00361 


Chr5 


- 


Zinc-ion binding 












1437432_a_at (L) 


Triml2a 


-186.81 


2.35E-07 


Chr7 


Idd27 


1443858 at (L) 


Triml2c (TrimS) 


— J 




unr/ 




1435665_at(L) 


Trim30d 


-3.73 


0.000726 


Chr7 


ldd27 


1421550_a_at(L) 


Trim34 


-2.46 


0.00119 


Chr7 


ldd27 


1451277_at(L) 


Zadh2 


-2.29 


0.00237 


Chris 


ldd2tl 


1434193_at(L) 


Zmym6 


-1.60 


0.000597 


Chr4 


IddU 


Protein/nucleic acid binding 










1424454_at (H) 


Tmeni87a 


9.56 


0.00354 


Chr2 


Idd13 


1455863 _at (L) 


Spata511 


-7.28 


0.000253 


Chr2 


- 


1444714_at(L) 


Dcdc2b 


-2.90 


0.00328 


Chr4 


- 


1431231_at(L) 


Histlh3f 


-2.33 


0.0028 


Chrl3 


- 


1435312_at(L) 


Paqr7 


-2.08 


0.00243 


Chr4 


- 


Otiier enzymes 












1442424_at (L) 


Ctdspl2 


-36.64 


0.000789 


Chr2 


Idd13 


1453009 _at (L) 


Cpm 


-12.24 


6.21E-06 


ChrlO 


- 


1418035 a at (L) 


Prim2 


-7.69 


4.78E-05 


Chrl 


Idd26 


1423569_at (L) 


Catm 


-5.45 


0.00467 


Chr2 


Iddl3 


1416494_at (L) 


Ndufs5 


-5.42 


9.99E-06 


Chr4 


- 


1427302_at (H) 


Enpp3 


3.32 


0.00172 


ChrlO 


- 


1427943 at (H) 


Acyp2 


2.53 


0.000285 


Chrll 


- 


144246S_a_at(L) 


Ppip5kl 


-2.84 


0.000837 


Chr2 


Iddl3 


1455075_at(L) 


Pigv 


-2.56 


0.00188 


Chr4 




1451998_at(L) 


Taspl 


-2.53 


0.0044 


Chr2 


Iddl3 


1433977_at(L) 


Hs3st3bl 


-2.48 


0.00247 


Chrll 




Otiier (unknown transcripts are indicated below in the legend) 








1444741_at(L) 


Dock2 


-31.19 


0.000307 


Chrll 




1436061 _at(L) 


Chafla 


-13.09 


0.000734 


Chrl 7 




1452426_x_at (L) 


LOC433762 


-13.07 


1.08E-05 


Chr4 




1452359 at (L) 


Rein 


-2.50 


0.00113 


Chr5 




1442824_at (L) 


Raly 


-12.31 


0.000974 


Chr2 


Idd13 


1442494_at(L) 


Ubr2 


-11.47 


6.75E-05 


Chrl 7 




1456635_at(L) 


SpllO 


-8.90 


0.00437 


Chrl 


Idd5.4a/t 


1456685_at(L) 


Nsg2 


-2.96 


0.0021 


Chrll 





Statistical analysis (1-way-ANOVA; p < 0.005, Benjamini-Hochberg) followed by hierarchical clustering identified 65 differentially expressed probe sets (representing 60 different 
genes) at 4 weeks; L and H, respectively, indicate lower or higher expression in NOD mice compared to control mice, NOR and C57BL/6. Fold change (FC) was calculated by ratio 
of means of expression in NOD mice versus controls. Dashes indicate those genes are not located within a known TID susceptibility region (Md); all indicated Mds (except one, 
ldd5.4a/5.4) were identified as conferring resistance to diabetes ( http : / /www.tl dbase.org ). Two unknown transcripts included: BE373131 (L; Chrl9) and A230098N10Rik (H; 
Chr5). Genes highlighted in bold font were differentially expressed at all 3 ages, 2, 3 and 4 weeks. 
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Table 5 

Enriched Gene ontology categories for tlie list of CD4 T-cell NOD altered genes. 



GO category 



2 wk (134 genes) 



3 wk (252 genes) 



No. of genes in 
category 



p-value 



No. of genes in 
category 



p-value 



4wk(185 genes) 



No. of genes in 
category 



p-value 



53 



52 
28 



Biological process 

Cellular 

metabolism 

Primary metabolism 

Nucleobase, 

nucleoside, 

nucleotide and 

nucleic acid 

metabolism 

Biopolymer 

modification 

Localization 

Cytoskeletal 

organization and 

biogenesis 

Transport 

Intracellular 

transport 

Secretion 

Leukocyte activation 
T cell activation 
Regulation of 
hydrolase activity 



Molecular function 

Catalytic activity - 
Hydrolase activity 22 
Hydrolase activity, 10 
acting on ester 
bonds 

Phosphoric ester 7 
hydrolase activity 

Hydrolase activity. - 
acting on acid 
anhydrides 

Pyrophosphatase - 
activity 

GTPase activity - 

Nucleotidyltransferase 
activity 

Inositol or phos- - 
phatidylinositol 
phosphatase activity 
Binding - 
Metal ion binding - 
Protein binding - 
Cation binding 28 
Transition metal ion 26 
binding 

Zinc ion binding 24 
Nucleic acid binding 28 



4.24e-5 



2.82e-5 
2.02e-3 



2.84e-3 
3.01 e-3 



1.32e-3 



8.09e-3 
2.49e-3 



6.04e-5 
4.70e-3 



8S 

86 



25 



43 
11 



36 
16 



40 
15 



10 

16 

15 

6 
5 



140 

55 

69 



1.13e-3 

7.06e-4 



6.53e-3 



1.13e-3 
9.07e-3 



7.49e-3 
1.64e-3 

5.70e-3 
8.80e-3 
3.17e-3 
8.07e-3 



7.88e-4 
6.23e-3 



1.99e-3 

4.48e-4 

1.12e-3 

3.49e-3 
7.76e-3 



3.87e-5 
2.85e-3 
1.89e-3 



67 



35 



1.88e-3 



8.88e-3 



60 
30 
14 



10 



2.89e-3 
2.30e-3 
S.65e-4 



1.75e-4 



99 
43 



39 
34 



29 



2.82e-4 



1.63e-4 



4.64e-3 
2.02e-3 



6.74e-3 
4.44e-4 



5.59e-4 



Cellular 

component 

Intracellular 

Intracellular 

organelle 

Intracellular 

membrane-bound 

organelle 

Nucleus 

Cytoplasm 

Endoplasmic 

reticulum 



65 

53 



48 



33 



2.81 e-5 

1.22e-3 



1.39e-3 



i.06e-3 



113 



80 



58 
15 



2.35e-7 

9.33e-4 



6.71e-4 



1.38e-4 
1.72e-3 



79 



5.64e-4 



One-way ANOVA (p < 0.05, Benjamini-Hochberg) of mRNA expression data followed by hierarchical clustering identified 134, 252, and 185 different genes at 2, 3, and 4 weeks, 
respectively, whose expression was altered in CD4 T-cells from NOD mice compared to two control strains (NOR and C57BL/6). These lists of NOD altered genes were analyzed 
in WebGestalt [http:/ /bioinfo.vanderbilt.edu/webgestalt) using the hypergeometric test (p < 0.01). Categories represented by > 4 genes are shown. Dashes indicate that the 
corresponding categories were not significantly enriched at the respective ages. Categories indicated in bold font were common to all 3 age groups. 
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Table 6 

Enriched KEGG pathways for the Hst of CD4 T-cell NOD altered genes. 



KEGG pathway 



Adjusted p-value 



No. of genes 



Genes in the pathway 



2 wk 

Metabolic patbways 

Purine metaboHsm 
Pantothenate and CoA biosynthesis 
Riboflavin metabohsm 
Nicotinate and nicotinamide 
metabolism 

Starch and sucrose metabolism 
Pyrimidine metabolism 

3 wk 

Metabolic pathways 

Riboflavin metabolism 

Purine metabolism 

Pantothenate and CoA biosynthesis 

Pyrimidine metabolism 

Nicotinate and nicotinamide 

metabolism 

Spliceosome 

Starch and sucrose metabolism 



l.Ole-5 

5.22e-5 
0.0001 
0.0001 
0.0003 

0.0008 
0.0039 



0.0005 

0.0008 
0,0008 
0,0008 
0,0013 
0,0008 

0,0029 
0,0031 



10 

2 
4 
2 
3 
2 

3 
2 



EnppS. NdufsS. CalntlO. Prim2, Catm, 

Dnmt3b, Akrlel, Enppl, Rrml 
EnppS, Prim2, Enppl, Rrml 
EnppS, Enppl 
EnppS, Enppl 
EnppS, Enppl 

EnppS, Enppl 
Prim2, Rrml 



EnppS, NdufsS, GalntlO, Prim2, Catm, 

Inppll, Dnmt3b, Akrlel, Dut, Enppl 

EnppS, Enppl 

EnppS, Prim2, Enppl 

EnppS, Enppl 

Prim2, Pnptl.Dut 

EnppS, Enppl 

Crnkll,Srnp200, Snwl 
EnppS, Enppl 



4 wk 

Metabolic pathways 

Pyruvate metabolism 



0,0002 
0,0008 



EnppS, NdufsS, GalntlO, Prim2, Catm, 

Pigv, Mdhl 
Acyp2, Mdhl 



One-way ANOVA (p < 0.005, Benjamini-Hochberg) of mRNA expression data followed by hierarchical clustering identified 58, 115, and 65 probe sets at 2, 3 and 4 weeks, respectively, 
whose expression was altered in CD4 T-cells from NOD mice compared to two control strains (NOR and C57BL/6). The lists of differentially expressed genes were analyzed in 
WebGestalt ( http://bioinfo.vanderbilt.edu/webgestalt ) for enriched KEGG pathways using the hypergeometric test (p < 0.01 , Benjamini-Hochberg). Pathways represented by >2 
genes are shown. Pathways or genes indicated in bold font were common to the 3 age groups. 



Table 7 

Transcription factor binding sites enriched in the promoters of CD4 T-cell NOD altered genes. 



Transcription factor binding sites 



p-value 



No. of genes 



Gene names 



2 wk 

Ar 



Irfl 



0.008 
0.015 



9 
11 



EnppS, TrimS4, Nrdl, Relll, Hisppd2a, 
Rrml, Spl40, Trib3, A530032D15Rik 
GalntlO, Trim S4, Khdrbsl, Pldn, 
Tmem87a, TrpSSbpl, Enppl, Mfap3, 
Psmb2, Rtn4, Snrnp200 



S wk 
Irfl 



lrf7 



ISRE 



3.74E-04 



0,002 



4.48E-04 



21 



17 



17 



GalntlO, TrimS4, Khdrbsl, Pldn, 
Tmem87a, TrpSSbpl, Trim21, Ahsal, 
Capn3, Enppl, lfi2711, Ifitl, Mfap3, Pi4k2b, 
Psmb2, PtpnlS, Rtn4, Snrp200, TnrcOa. 
1 700029101 Rik, 2510039O18Rik 
Khdrbsl, Tmem87a, Zfp69, Trim21, 
Trim30, Anxa6, Hisppd2a, lfi2711, Ifitl, 
Kiflb, Nsg2, Pdrgl, Pi4k2b, Pnptl, Zadh2, 
A1451 61 7, 4930402H24Rik 
EnppS, NdufsS, Triml2, Trim34, Khdrbsl, 
Zfp69,Trim21,Trim30, Dyxlcl, lfi2711, 
Ifitl, Pi4k2b. Pldn, Pnptl, Snrp200, Tnrc6a, 
A1451617 



4 wk 

Ar 



Irfl 



lrf7 



ISRE 



6.73E-04 



i.l3E-05 



2.50E-06 



5.58E-08 



11 



15 



15 



16 



EnppS, TrimS4, Nrdl, Relll, Hisppd2a, 
LyOcl, Ly86, Oasla, Spl40, Taspl, 
A530032D15Rik 

GalntlO, TrimS4, Khdrbsl, Pldn, 
Tmem87a, TrpSSbpl, Ifitl. Ifit3, Ly86, 
Oasla, Pigv, Rsad2, Rtp4, Stapl, Zfpl20 
Khdrbsl, Tmem87a, Zfp69, Trim30, 
Hisppd2a, Ifitl, lfit3, Kiflb, Nsg2, Oasla, 
Rsad2, Rtp4, Stapl, Zadh2, AI451617 
EnppS, NdufsS, Triml2. Trim34, Khdrbsl, 
Zfp69, Trim30, Ifitl, lfit3, Ly86. Oasla, 
Pldn, Rsad2. Rtp4, ZmymS, A1451617 



One-way ANOVA (p < 0.005, Benjamini-Hochberg) of mRNA expression data followed by hierarchical clustering identified 58, 115, and 65 probe sets at 2, 3 and 4 weeks, 
respectively, whose expression was altered in CD4 T-cells from NOD mice compared to two control strains (NOR and C57BL/6). The lists of differentially expressed genes were 
analyzed for enriched transcription factor binding sites in their promoters by using the PRIMA software of the EXPANDER suite (p < 0,05). Binding sites or genes indicated in bold 
font were common to the 3 age groups. ISRE - interferon sensitive response elements. 



40 



D.N. Kakoola et al. / Results in Immunology 4 (2014) 30-45 




Fig. 2. Molecular network generated by ingenuity patliway analysis (IPA) from the dataset of 2 week-old mice. The merged network was generated from the list of genes 
differentially expressed in CD4 T-cells from 2-week old NOD mice compared to both control strains (NOR and C57BL/6). The list was selected from the hierarchical cluster of 
362 genes that had highly significant (p < 0.005, Benjamini-Hochberg) expression differences between strains at 2 weeks of age. The genes derived from our uploaded gene list 
(known as focus genes) are represented by gray icons while genes (or endogenous chemicals) derived from the IPA knowledge base that could be algorithmically connected to the 
focus genes are represented by white icons. Different shapes of the symbols represent the different classes of genes/molecules, e.g. squares = cytokines/growth factors; ovals — 
transcription factors, etc. 



define novel TID susceptibility loci. These results attest to the utility 
of molecular network analyses in identifying novel Idds or genes not 
in known Idds that could act interactively upstream or downstream of 
Idd regions to contribute to diabetes development. The role of these 
genes needs to be confirmed in future studies. Several groups have 
determined that resistance genes in ldd9/ll (Chr4) regulate the di- 
abetogenic activity of CD4 T-cells [5,10,18,19]. Interestingly, two of 
the NOD altered genes (Khdrbsl and Ptp4a2) lie within WdS/H, with 
Khdrbsl (an adaptor protein involved in signal transduction cascades 
of several receptor systems, including T-cell signaling) being common 
to all 3 ages. In support of possible involvement of interaction of genes 
on several genetic regions in suppressing the diabetogenic activity of 
NOR CD4 T-cells, Chen et al. [5] reported that CD4 T-cells from NOR 
mice were somewhat more protective against diabetes than CD4 T- 
cells from NOD mice congenic for just the NOR-derived Idd9/ll. Our 
study provides support suggesting that resistance genes within lddl3 
(and their downstream genes) may act interactively with those in 
ldd9/ll (and possibly in unidentified Idd on Chrll) to regulate the 
diabetogenic activity of CD4 T-cells. 

In addition to the NOD CD4 T-cell altered genes discussed above, 
several other altered genes also lie within Idds (Tables 1-4). All (ex- 
cept Idd5.4a/5.4) have been identified as conferring resistance to 
diabetes ( http: / /www.tl dbase.org ). An interesting family of genes 
highlighted is the tripartite-motif (Trim) family. Compared to the 
leukocyte study, the current study expanded the list of Trim family 
members. These genes, whose expression was repressed in NOD mice, 
with that of Trims/ 12c and 12a virtually undetectable, all lie within 
Idd27 on Chr7. Trim proteins, which bear several domains, includ- 
ing three zinc-binding domains, constitute a family of ~60 molecules 



with diverse biological functions, including regulation of inflamma- 
tion and innate immunity [31,32]. The order of the domains is con- 
served throughout evolution supporting a common molecular prop- 
erty for these proteins. We propose that one (or several) of the Trim 
family members identified in our studies may play an important role 
in immune cells in initiation of autoimmune diabetes. To this end, 
Trim21, a gene that for a long time has been implicated in human 
autoimmune diseases Sjogren's syndrome and systemic lupus ery- 
thematosus, in which patients exhibit Trim21 autoreactivity [33], has 
been reported recently to regulate the innate immune response to 
intracellular dsDNA 134]. Another Trim (Trim28) has also been re- 
ported recently to be involved in the global regulation of CD4 T-cells 
[35]. Mice with a conditional T cell-specific deletion of Trim28 had a 
spontaneous autoimmune phenotype characterized by mononuclear 
tissue infiltration and altered cytokine balance leading to nonfunc- 
tional T regulatory cells (Tregs). 

Thus, the present study provides new insights into the genes/ 
pathways in CD4 T-cells that collectively might play a role in diabetes 
development - resistance in NOR mice or susceptibility in NOD mice. 
The BCL2L1 pathway and the genes within IddI3, Idd9/U and Idd27 
that are altered at all 3 ages (Pldn, Trp53bpl, Tmem87a, Ctdspl2, 
Gatm, Raly, Khdrbsl, and Triml2a and Trim5/12c) make particularly 
interesting candidate genes/pathways as they most likely represent 
basic genetic defects in the NOD mouse. While our study has high- 
lighted prime candidate genes, future studies will confirm which of 
the NOD altered genes and/or their interacting partners or regula- 
tors act interactively to effect the diabetogenic activity of NOD CD4 
T-cells. In our lab, we are investigating the genetic loci underlying the 
expression of the altered CD4 T-cell molecular network to define key 
regulatory loci. 

The majority of the CD4 T-cell NOD altered genes were repressed, 
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Fig. 3. Molecular network generated by ingenuity pathway analysis (IPA) from the dataset of 3 week-old mice. The merged network was generated from the list of genes 
differentially expressed in CD4 T-cells from 3-week old NOD mice compared to both control strains (NOR and C57BL/6). The list was selected from the hierarchical cluster of 
982 genes that had highly significant (p < 0.005, Benjamini-Hochberg) expression differences between strains at 3 weeks of age. The genes derived from our uploaded gene list 
(known as focus genes) are represented by gray icons while genes (or endogenous chemicals) derived from the IPA knowledge base that could be algorithmically connected to the 
focus genes are represented by white icons. Different shapes of the symbols represent the different classes of genes/molecules, e.g. squares = cytokines/growth factors; ovals = 
transcription factors, etc. 




Fig. 4. Molecular network generated by ingenuity pathway analysis (IPA) from the dataset of 4 week-old mice. The merged network was generated from the list of genes 
differentially expressed in CD4 T-cells from 4-week old NOD mice compared to both control strains (NOR and C57BL/6). The list was selected from the hierarchical cluster of 
581 genes that had highly significant (p < 0.005, Benjamini-Hochberg) expression differences between strains at 4 weeks of age. The genes derived from our uploaded gene list 
(known as focus genes) are represented by gray icons while genes (or endogenous chemicals) derived from the IPA knowledge base that could be algorithmically connected to the 
focus genes are represented by white icons. Different shapes of the symbols represent the different classes of genes/molecules, e.g. squares = cytokines/growth factors; ovals = 
transcription factors, etc.. 
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Table S 

Central genes in the Ingenuity Patliway Netwoiks of ttie CD4 T-cell NOD altered genes. 



Central Gene (IPA symbol) 



Gene symbol 



Total connections ( no. of Focus genes'' ) 



Connected Focus genes'" (shown as IPA 
symbols) 



2 wk 
HNF4A 



IFNG 

114 

TP53 

IL6 

MYC 

BCL2L1I' 

JUN 

IL15 

Prostaglandin E2 

NR3C1 

PARPl 

AFP 

EP300 

HRAS 



Hnf4a 



Ifng 

IL4 

Trp53 

1L6 

Myc 

BC12I1I' 

Jun 
1L15 

N/A 

Nr3cl 

Parpl 

App 

EpBOO 

Hras 



27(10) 



26 (2) 
26 (3) 
24(4) 
24(2) 
24(1) 
16(2) 
16(1) 
13(1) 
12(1) 
12(1) 
11(1) 
11(1) 
11(1) 
9(2) 



PLDN, NRDl, KHDRBSl, TMEM87A, 

SPATA5L1. RRMl. ALDH5A1, RBM39, 

SNAP23, C1ORF109 

CAPN3, Adaptor protein 1 

LCP2, GATM, Adaptor protein 1 

TP53BP1, PRIM2. RRMl, PLTP 

LCP2, PLTP 

DNMT3B 

GALNTIO, RTN4 

RBM39 

TRIMS 

ENPP3 

TP53BP1 

ASCC3L1 

HBA2 

SSI 8 

ENPP1,RTN4 



3wk 

INF 

TP53 

IFNG 

NFkB 
HNF4A 



TGFBl 
ERK 

Beta-estradiol 

p38MAPl< 

HA 

IGFl 

RELA 

IL15 

BCI^Ll" 

Prostaglandin E2 



Tnf 
Trp53 

Ifng 

N/A 
Hnf4a 



Tgfbl 

Ephb2 

N/A 

Mapkl4 

IL4 

Igfl 

Rela 

IL15 

BC12I1'' 

N/A 



30(6) 
30(7) 

28(5) 
28 (4) 
27(14) 



26 (4) 
25(5) 
25(6) 
23 (2) 
22(5) 
22(5) 
21(3) 
19(3) 
19(3) 
18(1) 



OASl. IFITIL, LBP, NC0A3, PDRGl, B1RC6 
TP53BP1, PRIM2, JDP2, PDRGl, BIRC6, 
DUT, NC0A3 

CAPN3, RTP4, TR1M21, 1F127, OASl 
ENPP1,LCP2, PNPTl.LBP 
PLDN, NRDl, KHDRBSl, TMEM87A, 
SPATA5L1, AHSAl, NC0A3, LBP, TMEM87B, 
FXCl, KIAAO409, CCDC21, PTPN18, PINKl 
PI4K2B, ENPPl, NC0A3, PlNKl 
LCP2, DNMT3B, ANXA6, INPPLl, LBP 
IFlTlL STC2, TCF15, NC0A3, PTP4A2, PlNKl 
LBP,JDP2 

LCP2, GATM, 1FIT1L,TRIM30, NC0A3 

KHDRBSl, INPPLl, NC0A3, TRPV2, MGA 

COMMD7,PNPTl,NCOA3 

TRIMS, FKBP1A,PTP4A2 

GALNTIO, RTN4, PNPTl 

ENPP3 



4 wk 



Interferon alpha 


Interferon alpha 


34(4) 


Oas, 1F1T1L,RSAD2,0AS1 


IFNG 


Ifng 


33 (6) 


RTP4, LY6C1, Oas, OASl, RSAD2, 1F1T3 


IL12 


1L12 


31 (3) 


TR1M30,LY6C1,1FIT1L 


114 


1L4 


30(5) 


LCP2, GATM, TR1M30, 1FIT3, IFITlL 


IL6 


1L6 


27(3) 


LCP2, 1F1T1L,LY86 


STAT4 


Stat4 


18(2) 


1F1T1LTR1M30 


111 5 


1L15 


17(1) 


TRIMS 


Prostaglandin E2 


N/A 


16(1) 


ENPP3 


TP53 


Trp53 


14(2) 


TPS3BP1, PRIM2 


BCHU" 


BC12I1I' 


14(1) 


GALNTIO 


HNF4A 


Hnf4a 


13(5) 


PLDN, NRDl, lOIDRBSl, TMEM87A, 
SPATASLl 


DDX58 


Ddx58 


11(2) 


1F1T3,1F1T1L 


SRC" 


Src'' 


10(2) 


KHDRBSl, PAQR7 


E2F4 


E2F4 


10(2) 


KHDRBSl, PRIM2 


IFIT2 


mt2 


9(1) 


1F1T3 



"Genes that came from the list of NOD altered genes. The connections to each gene (i.e. edges representing gene-gene relationships) in the molecular networks were manually 
counted. The central genes were then ranked from highest to lowest based on the number of total connections; the total # of connections represents the initial # of edges prior to 
editing of the networks for better image visibility. The table shows the 15 topmost central genes and their immediately connected focus genes. 

''BCL2L1 and SRC are located within the diabetes resistance region lddl3 (Chr2); BCL2L1 linked NOD altered genes (GALNTIO, RTN4 and PNPTl) are all located on Chrll (not in a 
known Idd) while SRC linked genes (PAQR7 and KHDRBSl ) are both located on Chr4 (KHDRBSl within the CD4 T-cell diabetogenic activity region ldd9/l 1 while PAQ.R7 is not in a 
known Idd). Central or focus genes/molecules indicated in bold font were common to all 3 age groups. N/A - not applicable. 



similar to the results of the spleen leukocyte study [25] and supported 
by many other studies in both mice and humans [36-42], Kodama et 
al. [36] reported global repression of genes in various tissues of NOD 
mice (including spleen cells) in comparison to NOD.BIO controls, a 
strain in which the NOD MHC haplotype is replaced with that from 
the nondiabetic BIO mice. Liston et al. [37] also found a global damp- 
ening in gene expression in NOD thymocytes of the genes involved in T 
cell negative selection. They identified several differentially expressed 



type 1 diabetes candidate genes, among them four genes (Ly6c, Prim2, 
Triml2, and Trim30) whose expression was also dampened in our 
study, thus providing supporting evidence for possible involvement 
of these candidate genes in CD4 T-cells as well. Heinig et al. [38], in 
a systems-genetics study investigating rat tissues and macrophages 
and human monocytes, discovered that the IRF7-driven inflamma- 
tory network (which is enriched for innate immune response genes) 
is associated with type 1 diabetes risk. Interestingly, homologs of 9 
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of the genes belonging to this network (Lcp2, Oasla, Rtp4, Ifitll, Ifitl, 
Ly6cl , Ifit3, Spll 0, and Trim21 ) were also repressed in NOD mice, sug- 
gesting a similar network may also be altered in CD4 T-cells of NOD 
mice and might contribute to diabetes development in the mouse. 
Similarly, in human TID studies, Elo et al. [39] found early suppres- 
sion of immune response gene expression in whole blood samples of 
children in the prediabetic phase who eventually developed clinical 
diabetes. Furthermore, Orban et al. [40] also found repression of all 
genes that were differentially expressed in whole untreated human 
peripheral blood CD4 T-cells from new onset TID patients, including 
genes involved in key immune functions, such as adhesion molecules 
lymphocyte function-associated antigen 1 (LFA-1) and P-selectin. Of 
note, L-selectin (the mouse homolog of P-selectin) was also of sig- 
nificantly lower expression in NOD compared to C57 (p = 3.2e-06 
and 0.0065 at 3 and 4 weeks, respectively). Functional analyses by 
Orban et al. further indicated that the CD4 T-cells were hyporespon- 
sive to stimulation and exhibited inhibition of cell cycle entry. In as 
much as the Orban et al. microarray study investigated whole un- 
treated CD4 T-cell populations rather than antigen specific suggests 
that the observed expression defect involves the polygenic CD4 T- 
cell population, and thus may signify a global CD4 T-cell repression. 
In support of this notion, other studies have also showed that CD4 
T-cells from TID patients have impaired activation (decreased pro- 
liferation) to non-specific primary but not diabetes-specific antigens 
[41 ]. It is worth noting that the human studies cited above used pe- 
ripheral blood samples providing further support for the potential 
usefulness of CD4 T-cells from spleen (a peripheral organ) in yielding 
results that may be of relevance to TID. Thus, future studies can test 
whether human peripheral blood cells (which are readily accessible) 
carry similar defects as those detected in our spleen CD4 T-cell study. 

Similar to the study by Orban et al., we also investigated whole 
untreated CD4 T-cells rather than antigen specific T cells. The fre- 
quency of T cells with a given antigen specificity is usually very low, 
estimated to be in the range of 1 or fewer in every 20,000 cells [40]. 
Therefore, we argue that the defect detected in our study also in- 
volves polygenic CD4 T-cells and is likely a basic genetic defect, in 
the least for those changes that were observed at all 3 age points. 
We cannot rule out the possibility that the differences obsereed, es- 
pecially the age-specific differences, may be due to influences other 
than genetic defects, e.g. activation of islet specific CD4 T-cells, es- 
pecially at 4 weeks of age. However, to this end, evaluation of our 
expression data revealed no evidence of an activated phenotype in 
NOD mice, suggesting that the changes in the gene expression may 
not be attributed to auto-reactive cells. For example, expression of 
two activation markers CD69 and CD38 was significantly lower in 
NOD (and NOR) compared to C57 (p = 8.26e-05 and 0.0056, respec- 
tively, at 4 weeks). CD69 expression was particularly interesting in 
that it showed no significant difference between the three species at 
2 weeks (p = 0.19) but became highly significantly upregulated at 3 
and 4 weeks in C57 (3-fold at each age, p = 4.57e-04 and 8.26e-05, 
respectively) compared to both NOD and NOR. Obviously this result 
was not reported in the "Section 3" as these types of changes were 
not the focus of this paper. Notwithstanding, it gives further support 
to the notion of a deficiency in T cell activation in the NOD mice, the 
effects of which NOR overcomes likely due to presence of resistance 
genes. 

Whereas the idea of a generalized repression of the immune sys- 
tem may seem counterintuitive (since later stages of TID are domi- 
nated by inappropriate activation of autoreactive T-cells), it may help 
explain previous observations that immunosuppression exacerbates 
autoimmune diabetes in NOD mice [42] while immunostimulation 
circumvents it [43]. Furthermore, it may not be surprising that rela- 
tive reduction of gene expression was seen at the early stages studied 
here, at which time tolerance would be induced. Of note, tolerance to 
self antigens requires the activation of self-reactive lymphocytes and 



their elimination by apoptosis as a result of this activation. Further- 
more, elimination of self-reactive lymphocytes in the periphery also 
requires activation of the regulatory mechanisms (such as regulatory 
T cells). Thus, both central and peripheral tolerances are active pro- 
cesses that require normal mitochondrial and metabolic function. Our 
gene Ontology and KEGG Pathway analyses (Tables 5 and 6), together 
with the leukocyte study [25], provided evidence for defects in mi- 
tochondria, metabolism, antigen processing/presentation and T cell 
activation/function and immune response. Furthermore, our prelim- 
inary functional studies investigating the mitochondrial/metabolic 
defect found impaired mitochondrial potential in NOD spleen leuko- 
cytes. All these data, together with the literature discussed above, 
support the idea of a global immune repression, which may lead 
to the breakdown of self-tolerance in autoimmune diabetes. Indeed, 
a recent study, using a mouse model of spontaneous autoimmune 
arthritis [44], suggested that efficient suppression of autoimmune 
diseases requires polyclonal regulatory T cell specificities rather than 
single antigen-specificities. Thus, we propose the following hypoth- 
esis. A genetic defect in metabolism/mitochondria results in a global 
repression of the immune system leading to a deficiency in immune 
tolerance, thus predisposing NOD mice to autoimmunity. Analysis of 
changes in gene expression and molecular pathways in NOD mice 
between different ages will shade further light on the defects that 
directly accompany initiation of insulitis, and subsequently develop- 
ment of diabetes. Furthermore, the defects in antigen presenting cells 
(such as B cells, macrophages and dendritic cells) may synergize with 
defects in the regulatory and effector T-cells to create dyshomeostasis 
in the early stages of autoimmune diabetes. Thus, investigation of the 
APC cell subsets is also warranted to provide a more comprehensive 
picture of the molecular pathophysiology of autoimmune diabetes. 

The promoter and molecular pathway analyses (Tables 7 and 8) 
identified several factors that may play a role in regulating the above 
discussed defect. Several of these factors (Hnf4a, Ifng, Trp53, Myc, 
111 5, Tnf, Tgfbl, beta-estradiol, 1L6 and Ar) were also identified by 
the spleen leukocyte study [25] indicating a strong involvement of 
the CD4 T-cells in the unfractionated immune cells. In addition, the 
current study identified novel regulators, including Interferon alpha/ 
beta, Irfl, Irf7, ISRE, III 2, and Stat4. Type 1 interferons regulate 3 
(Irfl, Irf7, ISRE) of the 4 transcription regulators identified in this 
study, suggesting a critical role for these cytokines in regulating gene 
expression abnormalities in NOD CD4 T-cells at the preinsulitis stage. 
Irfl is well known to control immune response gene expression 145] 
and has been demonstrated to be an essential element (in addition 
to Ifng and IL12) in the differentiation of naive T cells [46,47]. Irfl 
also functions as a transcription activator of the Tnf receptor and of 
genes induced by Ifng and type 1 interferons (including Irf7 and other 
ISGs) [45]. Together with the literature, our data provide support for 
a role for Irfl in regulating self-tolerance in autoimmune diabetes. 
Overall, our study captured new information, which, combined with 
future confirmatory studies, will facilitate a CD4 T-cell systems-based 
understanding of autoimmune diabetes and could ultimately lead to 
the development of novel therapeutic strategies. 
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